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We propose a method for determining the mass difference between two particles, h and h, that 
are nearly degenerate, with Am = m2 — mi <C mi. This method applies when (a) the li momentum 
can be measured, (b) h can only decay to h, and (c) li and h can be produced in the decays of 
a common mother particle. For small Am, h cannot be reconstructed directly, because its decay 
products are too soft to be detected. Despite this, we show that the existence of h can be established 
by observing the shift in the mother particle invariant-mass peak, when reconstructed from decays 
to h- We show that measuring this shift would allow us to extract Am. As an example, we study 
supersymmetric gauge-gravity hybrid_models in which Zi is a meta-stable charged slepton next-to- 
lightest supersymmetric particle and h is the next-to-lightest slepton with Am 5 GeV. 



Introduction. How well can we measure the masses 
of new particles at the LHC? This question might seem 
secondary at the early stages of discovery. It becomes sig- 
nificant, however, if there are two new particles that are 
nearly degenerate, with a mass difference of the order of a 
few GeV. In the busy LHC environment, it will be highly 
non-trivial to observe such mass differences. On the other 
hand, it is quite natural to expect such scenarios. If the 
new physics replicates the three-generation structure of 
the standard model, then we could easily have particles 
with the same gauge charges whose mass differences only 
depend on small "flavor effects." Obvious examples are 
the squarks and sleptons of supersymmetry, particularly 
the first- and second-generation ones. Here we propose 
an indirect method for resolving such mass differences, 
which applies when (a) the momentum of the lighter par- 
ticle can be measured, (b) the heavier particle can only 
decay to the lighter particle, and (c) the two particles can 
be produced in the decays of a common mother particle. 
For concreteness, we will discuss the case of two sleptons, 
h,2, with Am = TO2 ~- mi <C mi, with a meta-stable h. 
Such a spectrum is predicted by supersymmetric mod- 
els [H 0, that explain the masses and mixings of the 
standard model charged leptons and neutrinos in terms 
of broken flavor symmetries Q • More generally, this sce- 
nario is often realized in gauge-mediated supersymmetry 
breaking models [5*1 , and in large regions of the parame- 
ter space of models with gravity-mediated supersymme- 
try breaking 0]. Since the meta-stable li leaves a track 
in the muon detector, its momentum can be measured, 
and the event is fully reconstructible. Similar spectra and 
phenomena are also possible in other frameworks, such 
as for the Kaluza-Klein excitations of quarks and leptons 
in models with universal extra dimensions [7] . 

The method relies on the decays of a mother particle. 



in this case the neutralino XiJ to li^2- If we use direct 
decays to li to reconstruct the neutralino, the invariant 
mass distribution will be peaked at the correct neutralino 
mass. Some of the time, however, the neutralino decays 
to I2, which subsequently decays to h. The leptons pro- 
duced in the I2 h decay are relatively soft, with en- 
ergies typically of the order of Aw. They may therefore 
be lost, implying that the decays Xi h and Xi h 
have the same topology. Rather than blurring the pic- 
ture, however, we show that by attempting to reconstruct 
the Xi in both cases, one can in fact find two peaks: one 
at the neutralino mass M, and one slightly shifted to a 
lower value by an amount -Eshift ^ Am. Thus, measuring 
a shift in the neutralino mass peak will tell us that there 
are in fact two slepton states lighter than the neutralino, 
with a mass difference roughly given by i?shift- Further- 
more, if the two particles in question are scalars (which 
can be established from the angular distribution of the 
decays). Am can be determined in terms of i^shift, mi, 
and the neutralino mass M . 

The Shifted Peak. By assumption, the neutralino 
has two possible decays into sleptons. The first is the 
direct decay to h, 

X? - ifif ■ (1) 
The second is the decay to I2, 

X? - i^if , (2) 

followed by one of the two three-body decays [1, 0| 

l± ^ ifx^T , (3) 
It - IfX^^ , (4) 

where X^^ contains two opposite-sign (OS) leptons, and 
X^^ contains two same-sign (SS) leptons. Note that 
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the charge-flipping decays of Eq. Q resulting in SS lep- 
tons are possible because the neutralino is a Majorana 
fermion; SS leptons are also present in models other than 
supersymmetry when the decay is mediated by a vector 
boson or a scalar. We denote these lepton pairs by X 
to emphasize the fact that they are too soft to pass our 
cuts. Thus, the observed particles are the hard lepton 
from Eq. ([T]) or (l2|), and the long-lived slepton li from 
Eq. d), dSl), or dl. 

We can thus construct distributions for the following 
invariant masses-squared: 



"^1 = 



(5) 



where the li and I2 charges can be either opposite or the 
same. Obviously, 



(6) 



where M is the neutralino mass. However, because of the 
missing soft leptons. 



7^ M 



(7) 



and the peak of the mj^ distribution is shifted from M 

The identities of h and I2 depend, of course, on the fla- 
vor decompositions of h and I2, respectively. In one ex- 
treme case, if 1 1^2 are the left- and right-handed sleptons 
associated with the same flavor, the two leptons are iden- 
tical. Still, a peak in the invariant mass m?± ±, formed 

from events with SS sleptons and leptons, can only come 
from the decays of Eq. ([2|) and will therefore exhibit the 
shift i?shift- The analogous OS distribution will contain 
both types of events specified in Eqs. ([T|) and and so 
will exhibit a double peak structure, with the two peaks 
separated by -Eshift- 

In the opposite extreme, the two sleptons could be pure 
states of different flavors. In this case, the leptons li and 
I2 are different flavors, and there is no need to rely on the 
charges to separate the distributions. In any case, we can 
use the shifted peak to infer the existence of two distinct 
states. In the following, for the sake of simplicity, we will 
consider the second case, and take h = e and I2 = so 
that h = e and I2 = A*. 

We now turn to the calculation of the peak shift i^shift- 
We denote the four momentum of a particle a by Pa, and 
the h energy by Ei. Then, 



(8) 



From here on, we neglect the lepton masses. Working 
in the I2 rest frame, we take the x — y plane to be the 
plane of the muon and dilepton momenta, and px ■ We 
further take the muon direction to define the x-axis. The 



four-momenta of the hard muon and of the soft dilepton 
are then given by 



M' 



,t2 



2m2 



Px ^ "12 - Ei,-nJ El 



(9) 
(10) 



where n = {cos 6 , sin 9 , 0) . Substituting Eqs. © and (fTUl) 
into Eq. ([5]), we find 



Lfi 



TO2 



Ei-coseJEf-ml] , (11) 



with El and cos varying independently in the intervals 



El e 



mi, 



TO2 + rrii 



2m2 



cos 6* e [-!,+!] 



(12) 
(13) 



Let us now consider these quantities for Am <C toi. 
Working to leading order in Ato, we find that the maxi- 
mum value of El , given in Eq. (|12p , is 



TO2 -|- TO^ 
2TO2 



where 



mi 1 



Am 

TOi 



1 



We can therefore parameterize 



El = mi 1 



1 



(14) 



(15) 



(16) 



where < a < 1 varies from event to event. Note that 
the 0{x) piece vanishes. To leading order in the mass 
splitting, we then find 

m?^ - [(M^ + mj) + cosOy^iM'^ - ml)] x . 

(17) 

In reality, i?shift is defined by the peak of the m?^ distri- 
bution to be 



EsMit = M - ,Jm?^|pcak 



(18) 



and so it depends on the matrix elements governing the 
decays. Still, recalling that a < 1, the second term on 
the right-hand side of Eq. (fT7|) is always smaller than the 
first term by at least 



M2 



(19) 



Thus, a rough estimate for the mass splitting Am can be 
obtained from 



E, 



hift 



2Mmi 



-Am 



(20) 
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Since the ratio of Eq. (|T9|) will be measured, we will know 
the accuracy of this estimate. 

Beyond Leading Order. In fact, for a scalar I2, we 
can carry out the analysis exactly (including terms higher 
order in Am), and for all possible matrix elements, since 
we can determine the m-^ peak position based solely on 
kinematics. 

The crucial point is that, since the sleptons are scalars, 
the I2 — > hX decays are uniformly distributed in cosO. 
For a fixed Ei, then, these distributions are centered at 



2 2 



AI^ - mi 
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with width 



A(m? ) = 2- 



El 



(21) 



(22) 



m2 

As El decreases from {rn\ + rn\) / {2m2) to mi, the width 
gets smaller and smaller, until the distribution becomes 
infinitely thin and centered at 



m? = m 



y A'P - ml y KP + mim2 
' mi M 



Ar 



m2 



m,2 



(23) 

The total distribution is obtained by adding up all the 
contributions from different values of Ei with the appro- 
priate weights. 

We reach the following two conclusions. First, the peak 
of the total distribution is at m? . To see this, note that 



m? e 



m\ 



m2 
AP - ml 
m2 



El 



El 



El -ml,, 



El - ml 



(24) 



for all possible Ei. Thus every Ei contributes to this 
value of rri^^, and this is the only value of m?^ that every 
El contributes to. It is therefore the peak. 

Second, since the distribution is not symmetric about 
the peak, the mean of the distribution need not be at fn^^. 
This is relevant, because it means that the peak after 
experimental smearing need not be at However, the 
mean must satisfy 



ml e 



^2 2 , {AP-mDiml + mD 



(25) 



Thus we see that, to first order in Am/m2, the peak and 
the mean coincide, reproducing the result of Eq. (|^ . 

The Analysis. We now illustrate our method by 
simulating events in a concrete example model. Rather 
than simulate a realistic supersymmetric model, we use 
Herwig [13, Hi to specify a model to isolate the pro- 
cesses we are interested in, namely superpartners pro- 
duced through qq, qg, and gg pair production, followed 



by the cascade decays {g ^)q 



li^2- To do this. 



we choose mg = 650 GeV, all squarks degenerate with 



mg = 450 GeV, 



M = m^o 225.2 GeV 



m2 



mr = 139.9 GeV 

'2 



(26) 



mi = mr = 134.9 GeV 



and all other superpartners very heavy, so that they are 
decoupled from collider events. As mentioned above, it 
suffices for our purposes to consider flavor-diagonal soft 
terms and neglect left-right slepton mixing. We then let 
[i = CR and I2 = fj.R, and set B{xi h) = B{xi ^ ^2), 
as appropriate for the case of a gaugino x\ ^^'^ right- 
handed sleptons. 

We generate events for a ^/s = 14 GeV LHC and pass 
these events thr oug h a generic LHC detector simulation, 
AcerDET 1.0 df. We configure AcerDET as follows: 
electrons and muons are selected with > 6 GeV and 
1 77 1 < 2.5. Electrons and muons are considered to be iso- 
lated if they lie at a distance greater than Ai? > 0.4 
from other leptons or jets and if less than 10 GeV of 
energy is deposited in a cone of Ai? = 0.2, where 
Ai? = (A?])^ -f (A(/))^. The lepton momentum reso- 
lutions we use are parameterized from the results of Full 
Simulation of the ATLAS detector [l3|. (Our electrons 
are smeared according to a pseudo-rapidity-dependent 
parameterization, whilst muons are smeared according 
to the results for jryj < 1.1). AcerDET does not take 
into account lepton reconstruction efficiencies. We there- 
fore apply by hand a reconstruction efficiency of 90% to 
the muons and a reconstruction efficiency of 77% to the 
electrons. This gives 0.86 as the ratio of electron to muon 
reconstruction efficiency. We generate 20,000 events (be- 
fore any cuts or requirements are imposed), which cor- 
responds to 40,000 Xi decays. The supersymmetry cross 
section for the events of interest and for the parameters 
we have chosen is 50 pb, and so our event samples 
assume an integrated luminosity of £ ~ 0.4 fb^^. In a 
more realistic model, only a fraction e of all supersym- 
metry events will satisfy our event criteria, and so the 
assumed luminosity is 0.4 fb~^/e. 

The ATLAS trigger and reconstruction have been mod- 
ified recently to include the possibility of triggering on 



and reconstructing meta-stable sleptons |1J]. We con- 
servatively restrict ourselves to consideration of sleptons 
traveling fast enough to arrive in the principal time bin 
(25 ns wide). Furthermore, as mass resolution degrades 
as [3 approaches unity, it is necessary to place an upper 
limit on (3. Combining these two requirements we de- 
mand that the li candidates have velocity 0.6 < (3 < 0.8 
in agreement with Ref. [l^. Note that the upper cut on 
(3 will also greatly reduce backgrounds [l^l; we assume 
here that the remaining background events do not impact 
our results significantly. 

For the ^i, we take their momenta {pj^ \ from truth and 
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FIG. 1: Reconstructed h with 0.6 < jS < 0.8 after smearing h 
momenta with a Gaussian with cr = 0.05|pj-^ | centered on the 
true momentum and after smearing /3 with a Gaussian with 
a = 0.02. 



smear them by a Gaussian with a = 0.05 |p|J. We also 
smear the /? by a Gaussian with a — 0.02. These choices 
are again motivated by the results of Rcf. [15"| and their 
reconstruction of the li by measuring slepton time of 
flight with the ATLAS RPG chambers. We assume that 
mi is well measured, so we scale the four-momentum 
components to give the exact value of the mass. The re- 
sulting li mass distribution is given in Fig. [TJ We now 
take our measurement of mi to be 134.9 GeV. (Although 
our Gaussian fit actually gives 135.0 ±0.1 GeV, we as- 
sume that there is no systematic which would prevent us 
from obtaining 134.9 GeV exactly.) 

To eliminate the soft leptons from I2 decays, we impose 
a hard pT cut of 30 GeV. The px distribution of these 
soft leptons is shown in Fig. [21 from which we learn that 
such a cut is indeed a reasonable choice. Of course, in 
reality, the information presented in this Figure will not 
be available, and the experimenters will have to work 
by trial and error to find an optimal cut that gives the 
largest number of peak events while maintaining a clean 
peak. 

We first consider all OS l^e^ and l^^i"^ events and 
reconstruct the invariant masses mj^^ and "t-/^^, imposing 
the 30 GeV pt cut on leptons. Based on Eq. (PD)) . we 
expect mr 



mj , = 5.6 GeV: 



lie 

fhj 

hp. 



225.2 GeV 
219.6 GeV 



(27) 



The invariant mass distributions are shown in Figs. [3] 
and m We decompose each of these distributions into 
two pieces by fitting them to the sum of an exponen- 
tially falling contribution and a Gaussian distribution, 
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FIG. 2: The pr distribution of leptons produced by the three- 
body decays of h. 



with form 
dN 



dm 



il-fsi,)\a\e^^ + U 



2 1 

e 

TT cr 



(28) 

where a and mean have units of GeV ^ and GeV, respec- 
tively. These decompositions are also shown in Figs. [3] 
andlH The peaks of the Gaussian components then give 
us 



m 



ip 



225.4 ±0.1 GeV 
219.2 ±0.3 GeV 



(29) 



Thus, the separation between the peaks is experimentally 
well established. Using Eq. (|20p . we infer from the results 
of Eq. ((29)) that the selectron and the smuon are split in 
mass, with 



Am = 5.5 ±0.3 GeV 



(30) 



to be compared with our input value of 5.0 GeV. The 
exponentially-falling background in Figs. [3] and |4] is 
purely combinatoric. It consists of combinations of a li 
with a lepton on the "other side" of the event. In this 
toy MC we neglected other sources of SUSY background, 
from events with more than one lepton on each side of 
the decay. One can check however, (see Ref. [l3]) that 
these hardly affect our results, since the decay chain we 
considered here is the dominant one. 

So far we considered events with OS li and lepton. It 
is also instructive to examine the equivalent SS invariant 
mass distributions, shown in Figs. [5] and [6l As expected, 
we see no peak in the l^e^ distribution, and only the 
shifted peak in the l^^.^ distribution. Additionally we 
see that the muon peaks in the OS and SS samples are of 
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FIG. 3: The ife"^ invariant mass distribution. The fit param- 
eters a, a, /sig, and mean are defined in Eq. (|28p . 
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FIG. 5: The l^e^ invariant mass distribution. The fit param- 
eters a, a, /sig, and mean are defined in Eq. 
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FIG. 4: The iffi^ invariant mass distribution. The fit pa- 
rameters a, a, /sig, and mean are defined in Eq. (|28|) . 
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FIG. 6: The I'^/J.^ invariant mass distribution. The fit pa- 
rameters a, a, /sig, and mean are defined in Eq. (|28p . 



similar size. This implies that the probability for charge- 
preserving and charge-flipping decays are similar. This 
is not a given Q , and so provides an interesting addi- 
tional constraint on supersymmetric models. Note that, 
if the light sleptons are mixed flavor states, the SS dis- 
tribution can serve to reduce some of the background to 
the shifted peak. With mixed states, Fig. 0] would con- 
tain contamination from direct neutralino decays to li 
and a muon, but these must have opposite charges, and 
so do not contribute to the SS distributions. 

It would be nice to close this section by reporting a 
minimum Am that, given our assumptions about the ex- 



periment, could be established and measured. The power 
to discriminate peaks will depend on many things, how- 
ever. These include, but are not limited to, the cross sec- 
tion in each peak, the degree to which the peak shapes are 
modified by variation in acceptance across each plot, the 
cross section of the underlying backgrounds, and the un- 
certainties in those backgrounds. As mentioned above, 
with the strict upper cut on f3 that we used, we do 
not expect significant SM background. There would be 
however SUSY background as well as cavern back-splash 
background (not modeled in this investigation). For the 
numbers of events simulated in this study, and with not 
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all sources of backgrounds simulated, it is clear that the 
errors on the means fitted to the invariant mass peaks, 
being of order 0.3 GeV, would suggest that peak discrim- 
ination at the 5 sigma level is unlikely to be below 1 or 2 
GeV. Conversely, one would hope that discrimination at 
the level of 5 to 10 GeV (the width of the reconstructed 
invariant mass peaks) ought to be possible, and is indeed 
demonstrated here. 

Conclusions. We described a method for deter- 
mining whether a newly discovered long-lived particle is 
accompanied by another particle of almost degenerate 
mass, and for measuring the mass difference between the 
two. 

In the case of supersymmetry, which is the focus of our 
discussion, such measurements are central to analyzing 
the flavor structure of the theory, which will tell us about 
the origin of supersymmetry breaking. If the mechanism 
that mediates supersymm etry breaking is minimally fla- 
vor violating (MFV) [H, [ij, such as with pure gauge- 
mediation, then the mass splitting between the selectron 
and the smuon is expected to be determined by the muon- 
Yukawa squared, Am/m < ~ 10~^tan^/3, below the 
percent level. If we can experimentally establish a larger 
mass splitting, then we would obtain an intriguing clue 
for contributions that are not MFV fT]. 

If the dominant mechanism of supersymmetry break- 
ing is gauge-mediation, then the lightest charged slepton 
can be long-lived and leave a charged track in the muon 
detector of ATLAS/CMS. If the mass splitting is large 
enough that the transverse momenta of all the resultant 
leptonic decay products are frequently above the experi- 
mental reconstruction minimum (likely to be in the range 
5 to 10 GeV) then we can fully reconstruct the I2 h de- 
cays and measure the mass splitting directly. If, however, 
the mass splitting is only a few GeV, the leptons pro- 
duced in this decay are typically too soft to be detected. 
Despite this, as we have shown here, the hi invariant 
mass distribution will exhibit a double peak structure. 
The separation between the two peaks can be experi- 
mentally established and measured, and its value can be 
translated into a value of the slepton mass splitting. 

We note that a related study was presented in Ref. [IBl 
in the context of gravity-mediated supersymmetry break- 
ing, with a neutralino lightest supersymmetric particle. 
The events are characterized by missing energy, and 
the measurement of the slepton mass splitting is based 
on the kinematic edges of the distribution from 

With non-MFV supersymmetry breaking, we expect 
that the slepton mass eigenstates are not identical to the 
lepton flavor eigenstates. In this case, the decays we have 
considered contain a great deal of flavor mixing informa- 
tion i, [211,113, and our method can be used to measure 
not only the mass splitting, but also the mixing 17| . The 
SS sample will be particularly useful in this case, as it is 
only sensitive to neutralino decays to I2, and therefore 



only exhibits the "shifted peak." Counting the number 
of electrons and muons in this shifted peak gives a clean 
measurement of the flavor decomposition of l2- 
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